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An elementary approach to tangent space variation 
on Riemannian submanifolds * 

Jean-Daniel Boissonnat ^ Ramsay Dyer ^ Arijit Ghosh § 


Abstract 

We give asymptotically tight estimates of tangent space variation on Riemannian 
submanifolds of Euclidean space with respect to the local feature size of the submani¬ 
folds. We show that the result follows directly from structural properties of local feature 
size of the Riemannian submanifold and some elementary Euclidean geometry. We also 
show that using the tangent variation result one can prove a new structural property 
of local feature size function. This structural property is a generalization of a result of 
Giesen and Wagner [GW04i Lem. 7]. 

Keywords. Local feature size, Lipschitz function, Riemannian manifold and tangent 
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1 Introduction 

Let At be a compact Riemannian submanifold of the Euclidean space . For a point p 
in At, we denote by TpAt the tangent space of At at p. 

In this note we prove the following result: 

Theorem 1 (Tangent space variation) Let p, g G At and \\p — gll = tlfs(p). 

(i) If t <1/A then 

sinZ(rpAt, TqM) < tf{t), where f{t) = ^ ^ ' d) 

Observe that f{x) = § + 0{t), and if t < ^ then f{t) < 6. We will define Ifs(-) in 
Section n 

(ii) Better asymptotic constant: If t < then 

sinZ(rpAt, r,At) < 3t+ 0(th- (2) 

*A weaker version of this result first appeared in Chapter 4 of the third author’s PhD thesis [Ghol2l 
Lem. 4.3.2]. 
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A special case of this result in the context of 2-dimensional Riemannian submanifolds 
in was proved by Amenta and Dey [API 41 Thm. 2], and a weaker bound for the gen¬ 
eral case of Riemannian submanifolds of Euclidean space was proved by Niyogi, Smale and 
Weinberger [NSW081 Prop. 6.2 &: 6.3], using a global bound on the local feature size func¬ 
tion, rather than its local value. Belkin, Sun and Wang |BSW09l Lem. 3.4] showed that 
using cosine-law in the proof of the result of Niyogi, Smale and Weinberger one can show 
an optimal 0{t) variation bound. 

Our result is more general than these previous results, and we show that it is asymp¬ 
totically tight. Our proof is elementary and follows directly from structural properties of 
the local feature size function, Ifs(-), of Riemannian submanifolds and some elementary 
Euclidean geometry. 

We expect that the result of Niyogi et al. |NSW08l Prop. 6.2 & 6.3] together with Belkin 
et al.’s improvement |BSWn9[ Lem. 3.4] can be extended to the case of local feature size. 
Our main contribution is the simplicity of the proof of the tangent variation result wrt to 
the local feature size. We also show that using the tangent variation result one can prove a 
new structural property of local feature size, see Lemma [7| This result is a generalization 
of Lemma 7 from |GW04j . 

Notation. We denote the standard Euclidean distance between x, y € by ||x — y\\, 
and for a point x G and a set X C the distance between x and X will be denoted 
by d{x,X) = infj^ex ||x - y\\. 

For X G and r > 0, we denote balls and spheres by B{x, r) = {y G : \\x — y\\ < r}, 
B{x, r) = {y G : \\x — y|| < r}, and dB{p, r) = B{p, r) \ B{p, r). 

For S C we use B\g{p,r) to denote B{p,r) n S. We similarly define B\s{p,r) and 
dB\s{p,r). 

If U and V are vector subspaces of with dim([/) < dim(R), the angle between them 
is defined by 

Z{U,V) = max min Z{u,v) = arccos inf sup {u,v). 
u£U, v£V, u£U, vSiV, 

This is the largest principal angle between U and V. It is easy to show from the above 
definition that if dim(C/) = dim(P) then Z{U,V) = Z{V,U) (see |BG14[ Lem. 2.1]). 

The angle between affine subspaces is defined as the angle between the corresponding 
parallel vector subspaces. 

2 Medial axis and local feature size 

The medial axis of Ai is the closure of the set of points of that have more than one 
nearest neighbor on Ai. The local feature size oi x £ Ai, lfs(a:), is the distance of x to 
the medial axis of Ai |AB99) . As is well known and can be easily proved, Ifs is Lipschitz 
continuous, i.e., 

Ifs(x) < Ifs(y) -|- \\x — y||. 

Amenta and Dey [AD 141 Thm. 2], proved the following tangent variation result for the 
case of 2-dimensional Riemannian submanifolds of : 
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Theorem 2 (Two-dimensional case) Let M be 2-dimensional Riemannian submanifold 
o/M^. Letp, q £ Ai such that ||p —( 7 II = ilfs(p) with t < ^. Then sinZ(TpA^, TqM.) < 

The proof of the above result is restricted to the case where the dimension of the submanifold 
is two and the codimension is one. 

The infimum of Ifs over Ai is called the reach rch(Af) of Ai. Niyogi, Smale and 
Weinberger [NSWOSl Prop. 6.2 &: 6.3] first proved the following bound on the tangent 
variation on Riemannian submanifolds: 

Theorem 3 Let Ai be a Riemannian submanifold o/M'^. Letp, q £ Ai such that ||p—g] = 
trch(Wf) with t < 1/2. Then sinZ(TpWf, TqAi) < 2^yt{l — t). 

Belkin-Sun-Wang |BSW09[ Lem. 3.4] showed that using cosine-law in the proof of Theo¬ 
rem one can get the following improvement: 

Theorem 4 Let M be a Riemannian submanifold o/M^. Letp, q £ Ai such that ||p—(?! = 
trch(A4) with t < 1/2. Then sinZ(TpA4, TqAi) < 2tVl — = (2 — 0{t^))t. 

The proof of Theorem uses tools from differential geometry, such as parallel transport. 
Note that 0{t) bound in Theorem]^ is tight, see Section]^ But, as we have already men¬ 
tioned, the result use a global bound on the local feature size function, i.e. the reach, rather 
than its local value. We expect that this result can be extended to the case of local feature 
size function with more work connecting local feature size to other intrinsic properties of 
the submanifolds like the strong convexity radius. Our proof, on the other hand, is both 
elementary and works directly for local feature size function with no assumption on the 
global bound on the local feature size function. 

3 Proof of Theorem (i) 

The following lemma, proved in |GW04[ Lem. 6 & 7], states some basic properties of local 
feature size function. 

Lemma 5 1. Let p, q £ Ai such that \\p — gH = tlfs(p) with t < 1, then d{q,TpAi) < 

^lfs(p). 

2. Let p £ Ai and x £ TpAi such that \\p — x|| < tlfs{p) with t < \, then d{x,Ai) < 
2t^ Ifs(p). 

We prove Theorem using the above result. 

Let t = with t < 1/4. Using the fact that Ifs is 1-Lipschitz, we have 

(l-f)lfs(p) < Ifs(gr) < (1-I-f) Ifs(p) (3) 

Let u be an unit vector in TqAi. Let qu = q + tKs{q) ■ u, and let q'^^ denote the point 
closest to qu on Ai. Then, from Lemma I (2), we have \\qu — q'uW < 2t‘^lfs{q). Therefore, 
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Figure 1: In the figure a = Z([r, rii],TpM.) and A = d{ru,TpM). 


using the fact that Ifs(g) < (1 +1) lfs(j?) (see Eq. Q), we have 


\\p-qu\\ < \\p-q\\ + h-qu\\ + \\qu-(iu 

< tlfs(p) + (t + 2t^)lfs(g) 

< t Ifs(p) + (t + 2t^)(l + t)lfs(j?) 
= t (2 + 3t + 2t^) Ifs(p). 

Using Lemmas]^ (1) and (2), and Eq. Q, we have 


d{qu,TpM) 


< d{q'^,TpM) + \\qu - q'uW 
.2 

< - (2 + 3t + 2t'^) ^ Ifs(p) + 2tHis{q) 

f2 

< — (2 + 3t + 2t^) Ifs(p) + 2t^(l + t)lfs(j?) 
+2 

= — ((2 + 3t + 2t^)^ + 4(1 + Ifs(p) 


(4) 


Let r G TpM. be the point closest to q in TpAi, i.e., \\q — r|| = d{q,TpAi), and let v be 
an unit vector in TpM that makes the smallest angle with u. Let = r + tlfs(g) • u. Now 
observe that since ||g — r|| = \\qu — r„||, we have 


d{ru,TpM) < \\ru - qu\\ + d{qu,TpM) = d{q,TpM) + d{qu,TpM), (5) 

and the projection of the line segment [r, r^] onto TpM is parallel to v, which implies 
Z(n, v) = Z([r, ru],TpM). 
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Using the fact that d{r,ru) = tlis{q) > t(l —t)lfs(p) (from Eq. Q), d{q,TpM.) < ^lfs(p) 
(from Lemma[^(l)), and Eq. Q, we get 

smZ(u,t>) = sinZ([r, r„], TpA^) 

_ d{ru,TpM) 
d{r, ru) 

^ d{q, TpM) + d{qu, TpM) 

~ d{r, ru) 

~ 2d{r ru) + ^) + l) ^Hp) 

^ ^ f {2 + 3t + 2t‘^)‘^ + 4t + 5\ 

“V 2 - 2t J 

= tf{t) 

This completes the proof of Theorem (i). 

4 Proof of Theorem (ii) 

In this section we will give the proof of Theorem[^(ii) which improves on the the asymptotic 
constant in Theorem (i) from | to 3. We will use Theorem (i) itself to improve the 
asymptotic constant. 

Eirst observe that we can get an improvement on the bound given in Lemma 11(2): 
Lemma 6 Let p ^ M. and x G TpAi with \\p — x\\ = tlfs(p) with t < then 

d(x,M) < (1 - Vl -t2)lfs(p) = ^ 

By replacing Lemma[^(2) with Lemmaj^in the proof of Theorem[^(i) given in Section]^ 
we obtain the improved asymptotic constant stated in Theorem(ii). 

We will prove the following stronger result, which implies Lemma 

Lemma 7 Let iTp : A4 —>■ TpM denote the orthogonal projection of M onto TpM. Then 

(i) TTp restricted to B\j^{p,r) is an embedding, where r lfs(p)/10, 

(ii) B\TpM{P^'r') C 'Kp{B\_M{p,r)) where r' 191fs(p)/200, and 

(iii) let X G TpM with \\p — x|| = tlfs(p) with t < then ||x — 7r“^(x)|| < 

Proof Using the Morse-theory argument from |BC01l Prop. 12], we can show that 
B\M (Pt is an m-dimensional topological ball where m is the dimension of the manifold 
M. 
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(i) Observe that for all x and y G we have from Theoremand Lemma I (1) 


sinZ([x, y], TpM) < sinZ([x, y], T^M) + TpM) 

r 6r 

“ Ifs(x) Ifs(p) 

r 6r 

~ Ifs(p) — r Ifs(p) 

< 1 . 


since Ifs(x) > Ifs(p) — r 


( 6 ) 


The last inequality follows from the fact that r = This implies that iTp is a home- 

omorphism between r) and vTp (i3|_v((p, r)). Since if there exists x and y (/ x) in 

B{p,r) n such that 7rp(a:) = TTp{y) then Z([x, y],TpAi) = '7r/2 and that contradicts 


Eq. 


(ii) We now show that 

C 7rp{B\M{p,r)) where r'= (7) 

Since B^ji^{p,r) is a closed topolgical ball, and TTp restricted to B^j\^{p,r) is an embedding, 
it follows that TTp(B^_x^(p,r)) is also a topolgical ball containing p in its interior, and any 
point y G dTTp{B\j^{p,r)) is the image of a point 7r“^(?/) G dB\j^{p,r). Using Lemmaj^ (1) 
we find 

2 

\\p - y\\ > \\p - vr“^(i/)|| - ||7r“^(y) - y\\ > r - = r, 

and Equation Q follows. 


(iii) Let a: be a point in TpAf with \\p — x\\ = tlfs(p) with t < 2 ^. Since, from (ii), 
-S|Tpx(PA0 ^ '^p (-S|At(PA))) it follows that there exists x' G B\j^{r,p) such that 7rp(x') = 

x. Let ti = Erom Lemma (1) and the fact that vTp is a projection, we get 


||p — x'lP = tilfs(p)^ = Up — xll^ + ||x — x'lP < t^lfs(p)^ + ^lfs(p)^. 
Solving the above inequality we get 

f2 ^ _ f2 

< 1 _ = — i - , 

2 - i + VT^ 

and this, together with Lemma [^(1), implies 


d{x,M) < d{x,x') < ^lfs(p) < 


( 8 ) 


□ 
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Figure 2; Tangent variation on a sphere. 


5 Discussion 

Lower bound. It is easy to see that 0{t) bound on the tangent space variation on 
submanifolds is also tight. Consider a N — 1-dimensional unit sphere in For 

all p G Ifs(p) = 1. Let p, q ^ with \\p — g|| = t. Using elementary Euclidean 

geometry, see Fig[^ we can show that 



Using the above equation, and the fact that Z(TpS^ = Z{NpS^ ^), 

see [BG141 Lem. 2.1], we get the lower bound. 

Regarding the constants. The asymptotic constant we obtained in Theorem is 3, 
unlike the 2-dimensional case where it is 1, and we expect the asymptotic constant should 
be closer to 1. 
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